Collapsing Layers on Schwarzschild-Lemaitre Geodesics 
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We discuss Israel layers collapsing inward from rest at infinity along Schwarzschild- 
Lemaitre geodesics. The dynamics of the collapsing layer and its equation of state 
are developed. There is a general equation of state which is approximately polytropic 
in the limit of very low pressure. The equation of state establishes a new limit on 
the stress-density ratio. 
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I. INTRODUCTION 

Since their introduetion, Israel layers have played an increasingly important role 

in gravitational physics. Barrabes and Israel [3| began their paper with a description of the 
Israel layer as a thermodynamic phase boundary, but the initial applications of Israel layers 
considered metric matching in dynamic collapse processes involving dust shells, null shells, 

4j has summarized some of the early seminal 



and Hogan 



a. 



ing m dyi 



Poisson 



and cosmic string loops 

ff n 

work by Israel PJ,[2j, Barrabes [3|, de la Cruz [5[, Musgrave and Lake [6j, and Barrabes 



The Israel junetion conditions find wide application because they 

provide a simple dynamic boundary description for a variety of scenarios [8j. The use of 



layers as boundaries between Standard metrics has been summarized by Ansoldi [9j, and 



there is a growing literature that applies the thin she 



1 formalism to new areas such as 



shell quantum mechanics 9| - Jjj], layer dimensionality l^.[l^.jl5|, layers in extensions of 
general relativity [3] - [l9] and matching in perturbed spacetimes 20] - 24]. Astrophysical 
Problems using layers with constant spacetime character continue to be interesting, and 
include the investigation of collapse 25j, 26(, 27| and new phenomena such as gravastars 28] 
- 35j. A review of recent layer applications reveals increasingly complex layer models, the 
gravastar models, for example, involving multi-layer construetions, and there is an increasing 
interest in describing more realistic layers [36| . [371] . 
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A layer metric is determined by the metrics of the two bounding manifolds. The stress- 
energy of the layer is determined by jumps in the extrinsic curvatures of the metrics on either 
side of the layer. For example, the Israel layer between exterior vacuum Schwarzschild and 
interior Minkowski has a 2+1 metric as seen from both sides of the layer 4|] 

ds 2 = -dr 2 + R 2 (r)dü 2 . (1) 

The development of R{t) tracks the dynamics of the layer and, as in any general relativity 
Solution with stress-energy, the equation of state of the layer is important in characterizing 
the dynamics. 

In this work we discuss layers which start from rest at infinity in an exterior Schwarzschild 
metric with mass parameter m and drop inward along a geodesic. The Lemaitre form of the 
Schwarzschild metric is adapted to this particular geodesic motion but the layer motion can 
be simply described with the usual Schwarzschild metric. The interior space is Schwarzschild 
with mass parameter M. Although the motion of the layer is simple, the equation of state 
is cubic in the stress, quartic in the density and it provides a new restriction on the ränge 
of the stress/density ratio for collapsing layers. 

Layers with simple motions, such as the layer considered here, can require more physical 
content than the density and pressure coming from simple perfect fluid or polytropic mod- 
els. We quote from j^öj]: "... models satisfying a closed equation of state are relatively few 
and essentially restricted to dust or linear barotropic models. However, it is worth stressing 
the importance of considering general equations of state in realistic models, describing for 
instance, very high density regimes for white dwarfs and neutron stars." Layers with physi- 
cally realistic equations of state are important in studying the collapse scenarios that result 
in objects like neutron stars. 

In the next section we describe the time development of the radial function, R(t), in the 
exterior Schwarzschild metric. The stress-energy, the equation of state of the layer and its 
symmetries are discussed in Section III, and we close with a Discussion. Extrinsic curvatures 
are computed in an Appendix. 
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II. LAYER GEOMETRY 

A. Bounding metrics 

The 2+1 layer is bounded by two metrics 

ds 2 = g ab ±dx a dx b = -f±dt 2 + {l/f±)dr 2 + r 2 dü 2 (2) 
with +/- denoting exterior/interior. 

/+ = 1 - 2m /r (3a) 
/_ = 1 - 2M/r. (3b) 

The metrics correspond to an itlq Schwarschild exterior and an M = const Schwarzschild 
interior. The layer is described by r = R(t) and t = T±(t). The velocity of the layer as 
seen by observers moving with the layer in the bounding spaces is U± — (T±, R, 0, 0), with 
corresponding normal vector n a ± = (— R, T±, 0, 0) Q]. The normal vector is chosen to be 
outward pointing, toward the exterior spacetime. The velocity normalization imposes the 
condition 

MBfc-iim- 1 - (4) 

This normalization is used to dehne function ß± 

ß 2 ± :=fl{R)t 2 ± = f ± {R) + R 2 . (5) 

ß± is chosen positive with T± > and R < describing an infalling layer. An expanding 
layer can be described by proper time inversion. The equation of state of the layer will not 
depend on the direction of motion. The induced metric on the layer is 

with 

^ = _ ( «k) V± + (7) 

h TT ± = — 1 from the velocity normalization. The metric of the layer from both sides of the 
boundary is given by Eq.(TfJ). 

In the next section, R(t) will be fixed by the requirement of geodesic layer motion in 
the exterior spacetime. Since both sides of the layer will agree on the layer metric, this will 
determine the interior radial function. 
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B. Exterior geodesics 



For Schwarzschild geodesics in the exterior spacetime we have 

dR 



— = ±y/E% - 1 + 2m /R (8) 
dT+ _ E Q 

dr ~ 1- 2m /R 1 ' 

where Eq describes the initial point of the layer motion. Eq = 1, dR/dr < 1 corresponds 

to a layer beginning its inward drop from infinity. The relations to use in describing the 
infalling layer from the exterior spacetime are 

^ = -Ä/Ä (10) 



dT: 1 



dr l-2m /Ä" 



(11) 



Integrating, one finds 



i?(r) = (2m ) 1/3 [(3/2)( Cl -r)] 2 / 3 , (12) 
which completely determines the layer dynamics. The metric of the layer is 

ds 2 = -dr 2 + [(3/2)V2^] 4 / 3 ( Cl - r) 4 /W. (13) 

III. LAYER STRESS-ENERGY 

A. Fluid description 

We assume a perfect fluid stress-energy for the layer. The co-moving velocity is U % = 
(1,0,0) with density a and stress-energy 

:= (<r + P)U% + Ph'j. (14) 

The stress-energy of the 2+1 layer is related to jumps in the extrinsic curvatures of 
the bounding metrics 

-871-5^ := < K) > - < K > h) 

SnS; = -<K e e +K*> (15a) 

-8nS e e = - <K^ + K;> (15b) 

-8ttSJ = -<K e 9 +K; >, (15c) 
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where = [—l,R 2 ,R 2 sm 2 6], h % - = and where 

< K] > : K) . K) (16) 

is the jump in the extrinsic curvature and K := K\. The extrinsic curvatures are worked 
out in Appendix A. For the constant mass exterior Schwarzschild metric they are 

<- = 4+ = l / R ( 17 ) 

K + = o. (18) 

In the interior, with mass parameter M, the extrinsic curvatures are 

Kl = K*_ = ß_/R (19) 

M — mn , , 

= HPT- < 20) 



The layer density is 



where 



4™=-<A' e 9 >=^^ (21) 



/L^l+2^. (22) 

For i? > 0, the ränge of ß_ is restricted to /?_ > 1 for positive layer density. This requires 
mo > M. The stress is 

-ß- + 1 mo - M 



5ttP = < + KT > 



R R 2 ß- 



8,P = -^- 1 ) 2 i. (23) 



A layer mass m^ can be calculated from the layer density 

rriL = 4:iiaR 2 

The physical picture has an exterior Schwarzschild metric with mass parameter mo formed 
from the stress-energy of the layer plus the interior Schwarzschild metric with mass M. The 
binding energy of the layer contributes to the value of m , as can be seen by eliminating ß_ 
between equations fl22|) and ff24|) : 

2 

m L = m -M-^. (25) 
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The layer mass is not constant due to work done by the stress 



d(oR 2 
dr 



-P 



dR 2 
dr 



(26) 



This has also been noted by Carr and Yahil 39| for a self-similar Friedman universe and by 
Visser 40j in dynamic wormholes. 



B. Equation of State 

The relation between the density and stress for Lemaitre layers is 

p=-^r (27) 

or 

ß - ~ TTTpfa (28) 

Expression (1251) can be rewritten as 

ß _i = _ ip i" ß + 1 = 2 + 4J > 

or 



[1 + AP/a) 2 ' 

Eliminating R between Equations (|21ß and (1221) allows a to be expressed as 

°" 87r(m - M) 

Substituting for (/3_ — 1) and (ßt — 1) in terms of P/a above, yields the equation of state 
for the Lemaitre layer 

2 _ l n(m -M) 1 a 3 (l + AP/af 
P ~ [ 4 J (l + 2P/a) • (29) 
This expression can be rewritten for the density or stress in terms of the ratio x := P/a. 

4x 2 (l + 2x) 

7r(m - Af)(l + 4x) 3 ' 1 ' 

For a static membrane under tension, x is the negative of the sound speed squared, dPj da. 
The 'x' ranges of possible interest are 

I :0<x (31a) 

II : -l/4<x <0 (31b) 

Ii7 : -1/2 < x < -1/4 (31c) 

IV : -1 < x < -1/2 (31d) 
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The relative size of m to M provides three mass regions to consider. m = M has no 
extrinsic curvature jumps and no layer. For m > M, ß_ > 1 is required for positive layer 
density. This eliminates the positive ränge of x. ß- > requires — 1/4 < x and this eliminates 
the third and fourth region. Region II is the only allowed ränge and this sets a new limit 
on the stress/density ratio for layers collapsing in a Schwarzschild geometry. Figure 1 shows 
the Variation of the stress and density in Region II. A portion of Region I has been included 
for clarity. 
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FIG. 1: Scaled density and pressure vs pressure/density 




In general the equation of state does not have any of the usual simple forms. It is a cubic 
equation in pressure and, depending on the value of er, can contain simpler equations of state 
in the root structure. Rewriting Eq. (1291) . one finds 

p := p(7r/4)(m - M) 
ä:=a(vr/4)(m -M) 
JJ 3 (64ä - 2) + p 2 (48ä 2 - ä) + Yl^p + ö 4 = (32) 

The equation of state, Eq. (l32|) . is plotted in Figure 2 showing the regions where, for a Single 
density value, there are three real roots, two real roots, and one real root giving associated 
values of the pressure. The negative pressure region corresponds to the layer tension. In this 
region, for small density, there are two possible tensions for each value of density while for 
larger density, there is a 1-1 nondinear relationship. 
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FIG. 2: Scaled density vs scaled pressure for the equation of state 



Expanding Eq. (l29l) one finds 



P ~ —, - 10a; + ...) 

7r(m — M) ' 



For x « 1 the equation of state appears approximately polytropic with index 2 41]. 



Layers are most commonly described as perfect fluids or scalar fields. Within the perfect 



fluid models, several EOS 



251 ] , [38J , [42j have been suggested to analyze the layer stress energy 



content. A commonly discussed form is 

P = -na (33) 



which plays a crucial role in 3+1 homothetic Solutions [43j. The layer is dropping on a 
Schwarzschld/Lemaitre geodesic. The actual equation of state does not change but it can 
be modeled with simpler forms having time dependent parameters. This is a different point 
of view than many dynamic layer treatments which assume a static equation of state and 
let that drive the layer dynamics. For example, Eq. (1271) for the Lemaitre layer implies that 
n(r) changes as the layer drops with the n value given by 

n(r) = -x = ß -^. (34) 

n(r) runs through the ränge < n < 1/4 for the positive density layer. We note that x = 
is the dust Solution corresponding to m = M, the vanishing layer condition, or to very 
large R. Positive ß- > 1 describes an inward moving layer with positive layer density, mass 
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and tension. The constraint m > M implies that r = M is inside the r = 2m trapped 
surface. The physical picture of the layer dropping inward toward an interior object of mass 
M implies that the layer will cross the 2mo horizon. The value of x at the horizon, xh, is 
parameterized by the ratio M/mo 

4x H = - 1 - 1. (35) 

y/2 - M/m 

The maximum value of xh, %h ~ —0.07, occurs for a Minkowski, M = interior, with xh 
decreasing as the size of the interior M approaches rriQ. With M = 0, the tension of the 
layer is P = —.07a at R = 2m , increasing inward to —1/4 at the R — singularity. The 
layer mass at R = 2mo is less than mo, since mo will include the layer binding energy. 

The region, < M, describes a layer with negative density and mass and with ß_ < 1, 
with x positive. For this case R = 2mo lies inside R = IM. From Eq. (l22|) . ß_ = 
at R = 2(M — mo). At this radius, the layer has infinite tension. The interior Schwarzschild 
Solution could be the exterior of a third core fluid region with a boundary at R > 2M, 
eliminating trapped surfaces from the vacuum motion. 

Another indication of the complexity of the Lemaitre layer is the sound speed, 

da 

Using the equation of state, this is 



Vg is negative over the ränge — 1/4 < x < indicating instabilities in the layer 44]. 



IV. DISCUSSION 



In this paper we have discussed the geodesic collapse of an Israel layer in an exterior 
Schwarzschild metric. The geodesic drop of the layer determines the dynamic development of 
the layer's motion. Although their geodesic motion is simple, Lemaitre layers are physically 
interesting structures because of their complex equation of state, which can approximate 
some of the simpler equations of state commonly used in layer dynamics. Some of the 
simplest equations of state can be described by time dependent models of the actual stress 
density relation. A barotropic equation of state is often used to describe Israel layers, for 
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example in effective potential treaments [26j , [45] , [46j and discussion of bubble structures 



4=71 ] . In the simple model considered here, the layer can have a time dependent barotropic 
relation, P = —n{r)a. It has an R — > limiting value P = — er/4, with n(r) still very small 
as the layer reaches the R = 2mo surface. 

The ränge of P/a for positive density, has a lower bound of —1/4. This is a new restric- 
tion on the ränge of the stress/density ratio for collapsing layers. Garfinkle and Vuille 26] 
considered the layer between two Schwarzschild metrics with no motion restrictions. The 
value P/a = —1/2 was identified as a boundary for their effective potential 45j, 46] to have 
stationary points. The effective potential for the Lemaitre layer is simple with no stationary 
points or turning points and the boundary P = —er/4 occurs at the point R — 0. 

Fluids collapsing along geodesics in a Schwarzschild metric have been considered by 



Lynden-Bell and Lemos 



48j , generalizing the description of a Newtonian cold gas collapsing 



from rest at infinity 49]. The collapse of the Lemaitre layer and the cold gas have some 
similarities. They are both described by the the same R(t) and for either case Eq.( fl2l) can 
be written as 

4mp R^r) 3/a 
3 2m 

where Tj is the proper time to reach Ri and r c is the time to reach R = 0. In the less limiting 
case of a Newtonian sphere with inferior mass m, Lynden-Bell and Lemos show that r c can 
be written as a function of the mass 



r c = Bm + To, 

where Tq is a constant of integration. Generalizing this relation to Einstein's gravity, b = 1 
is identified as a value of special interest corresponding to metric homothetic symmetry. For 
the Lemaitre layer, r c is just the constant c\ = r , corresponding to B = in the cold gas 
Solution. r c is related to the time th for a layer to reach the R = 2m surface by 

t c = 4m /3 + th- 

The Lemaitre layer can be considered as a special member of the b — 1, B = geodesic 
Solutions which include cold gas collapse. One also notes that the exterior bounding Lemaitre 
metric exhibits second kind kinematic similarity 50(, reflected by R(t) in Eq. (fl2]) . which is 
transferred to the layer and maintained thoughout the motion. 
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There is growing interest in Israel layers as objects of research rather than an aid in 
matching metrics. The increasing interest in expanding the physical properties of Israel 



layers has led to descriptions incorporating layer parameters like shell thickness 36], 37] 



and surface packing fraction [131 ] . While Israel layers are still linked to the properties of 
the two bounding spaces, even simple layers such as the Lemaitre layer considered here, 
require more physical content than the density and pressure coming from simple perfect fluid 
or polytropic models. The time dependent models for the actual stress density relations 
described here, represent one way these simple models can still be used effectively and related 
to the literature using simple equations of state. 

APPENDIX A: EXTRINSIC CURVATURES 

The stress-energy of the layer is determined by jumps in the extrinsic curvatures, K a b of 
the bounding metrics across the Israel layer. The bounding metrics considered here are 

ds 2 = -f±dt 2 + (l/f±)dr 2 + r 2 dü 2 (AI) 

with (/+ , /_) = (1 — 2mo/r, 1 — 2M/r). The layer is described by r = R(r) and t = 
T±(t). The velocity of the layer as seen by observers moving with the layer in the bounding 
spaces is U± = (T±,R,Q,Ö), with corresponding outward pointing normal vector n a ± = 



(— R, T±, 0, 0) {4]. The velocity normalization imposes 

MBft-fSs)- 1 - <A2) 
The normalization is used to dehne a function, ß± 

ß 2 ± :=fl(R)f 2 ± = f ± (R) + R 2 . (A3) 

We choose ß± > 0. ß± < corresponds to r — ► — r. The metric of the layer from both sides 
of the boundary is 

ds 2 = -dr 2 + R 2 (r)dü 2 . 

The extrinsic curvature is defined in terms of the normal to the layer, n a and the projection 
Operator onto the layer, h a h as 

Kn := n a , b hth) (A4) 
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The two extrinsic curvatures are K e e and K T T . K e e is easily calculated from the definition and 
one finds 

Ät .^.^* (A5) 

where we use /3_ to indicate the ränge of ß± to be used. Using the infalling geodesic value 



for R = — sj 2m / R we have 



R R y ' 

K T T is more easily calculated from the acceleration K T T = n l a,i. The acceleration vector is 

a* = [f+Tj-, Ä+X, 0, 0]. (A7) 
Dötting with the normal and using the definition of ß±, one finds 

K± = % • (A8) 
All cases will have the same exterior functions. In the inferior, the extrinsic curvatures are 



ß ^ l+ 2^-M) (A9) 

Kl_ = ß_/R (A10) 
K.-^. (AU, 
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